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The first Townsend coefficient for Ar-C02 based gas mixtures has been measured over 
a wide range of reduced electric field. The experimental setup and the measurement 
technique are described here. A linear superposition model has also been successfully 
applied. 
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1. Introduction 

The Townsend coefficient, defined as the number of ions produced per unit path by 
a single electron traversing a gaseous medium, is of fundamental importance in all the 
discharge processes, in particular in the description of the electronic gain of gas ioniza- 
tion detectors[T] or plasma discharge devices^. Direct measurements of the Townsend 



coefficient are poorly reported in the literature |3|4|5|7|8] . therefore the most used esti- 
mate of the Townsend coefficient is based upon a numerical solution of the Boltzmann 
equation |9llUj . 

In this paper we report the measurement of the effective Townsend coefficient for Ar- 
CO2 gas mixtures in the typical working condition of high gain MWPC 0. We have 
derived the effective Townsend coefficient a from the measurements of the gas gain M in 
a cylindrical test tube, using a new method that is briefly described in the following. 

The gas gain M, defined as the ratio of the anode current / over the primary ionization 
current Jq that would be measured operating the tube in the ionization chamber mode, 
can be calculated integrating the Townsend coefficient from the starting point of the 
avalanche ro to the surface of the wire^I]. In practice we have 

InM = r a{S)dr (1) 

Jro 

where S is the reduced field, which in cylindrical geometry will be 

S{r) = . , (2) 

being V the voltage applied to the anode, the anode wire radius, the cathode radius 
and p the pressure of the gas. It is worth noticing that the coefficient a which appears in 
Eq. is an effective coefficient in the sense that it is the difference between the absolute 
coefficient and the attachment coefficient. Changing the integration variable [7] and using 
Eq. (j21), we can recast Eq. (Q) in the form 

where Sa = S{ra) and Sq = S^vq). A small change in V is equivalent to a small change 
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in the value of Sa- Differentiating Eq. Q respect to Sa, treated as an independent 
variable, we obtain after some algebraic manipulations 



which we have used to extract from our measurements the Townsend coefficient. It is 
noteworthy that, being always InM -C dlnM/dlnV, Eq. (@)) can be used to estimate a 
even if the primary ionization current Jq is poorly known, because o? In M = o? In / if Jq is 
constant. 

2. Experimental setup 

The experimental setup used in this work is shown in Fig^ The active counting volume 
is a cylinder with a thin anode wire of gold plated tungsten of 30 fim diameter, soft 
soldered to the bronze feed-throughs and subjected to a mechanical tension of 60 g. The 
wire is accurately located in the center of a precision stainless steel tube of (5.64 ± 0.01) 
mm inner diameter and (0.18 ± 0.01) mm wall thickness. Research grade pure gases were 
mixed, with relative percentages controlled to 0.5% accuracy level by two computer driven 
mass flowmeters; during operation the test tube has been fluxed with a total flow of 10 
seem. The high voltage is given by a power supply device, controlled with a IEEE-488 
interface, while the current is read with a Keithley Mod. 485 picoammeter, using the 
internal logarithmic conversion and remotely controlling the fullscale. The calibration 
error of the picoammeter stated by the constructor is 0.4% plus 1/2 LSD of the readout. 

The active volume of the tube has been exposed to a radioactive ^^Fe source, whose 
intensity was 100 /iCi. The X-rays reach the gas volume by passing through a small 
mylar window of 50 /xm thickness placed on the cylinder external surface. The distance 
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Figure 1. Schematic view of the experimental set up 

of the source from the mylar window has been adjusted to give a rate of current pulses 
from X-rays interactions in the active volume of the tube of about 100 kHz. The primary 
ionization current Jq has not been directly measured. Instead it has been calculated from 
the above quoted rate multiplied by the average primary charge released in the counter 
by the 5.9 keV X-rays, taking into account the escape for argon. As we said in the 
discussion of Eq (0)) the uncertainty in the determination of Jq has practically no effect 
on the calculated value of a. However we observe that also if the flux of X-rays from the 
radioactive source is constant, the counting rate will depend upon the relative composition 
of the gas mixture, due to the different absorption of X-rays in the two gases, and the 
absolute pressure and temperature which determine the density of the gas. In fact the 
primary ionization current Jq will be 

Jo oc (i-e-^*-^'-^-*' ) . 

Exil-Pesc) + iEx-E^-)P,sc ^ 
WAr 

+ f 1 - e-'^^oa ^^o^ (5) 

^ ^ WC02 

where /i are the respective mass absorption coefficients in cm^/g, for an X-ray with energy 
Ex, E^^ is the shell K ionization potential for argon, Pesc the escape probabihty of the 



fluorescence X-ray produced by radiative deexcitation of tlie K-sliell, and w tlie average 
energy for producing one ionization pair in tlie respective gas. Finally the thickness Xgas 
of the respective gas in g/ cm^ is given by the well known formula 



^,as = ^^^fl^(t) (6) 

being (t) the average tracklength of the X-ray in the gas, rugas the molecular weight of 
the gas, Pgas its partial pressure, T the absolute temperature and R the gas constant. We 
observe that a small variation of the pressure and/or temperature of the gas is amplified 
by the exponential dependance in Eq. (0). In order to correct for small changes of the 
pressure and temperature of the gas mixture during the run, we have monitored the 
pressure of the gas with an accuracy of ±0.5 mb and the temperature of the test tube 
with an accuracy of ±0.1 C. 

3. Experimental results 

We report here the results obtained in five different experimental runs, repeated in 
order to check the reproducibility of the measurements. Each run has been performed 
in completely automatic mode. We have adopted a rather complex strategy after some 
preliminary measurements, that have given to us a feeling for the various dynamical time 
scales involved in the system. First we have observed that the flow of the gas stabilizes 
at the assigned partial flows only after a period of at least half a day. During this period 
we found that a constant applied voltage of 1000 V helps the conditioning of the test 
tube to the new mixture. After this time the anode voltage is changed by the program 
to the minimum voltage and raised in steps of 25 V, with slow ramp of 2 V/s. When the 
high voltage has reached the programmed value, the program allows for additional 180 s 
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Figure 2. Townsend coefficient vs. reduced electric field. The six curves are obtained for 
Ar-C02 mixtures spanning from 100-0% to 0-100% relative percentages. The thin dot- 
dashed lines are the predictions of Eq. (jS)) with the parameters A and B have been fitted 
separately for the pure gases. The heavy dashed lines are obtained from a two-dimensional 
fit of the entire data set with the same formula (see text). 

waiting time, before starting the current measurements at this voltage. The measurements 
are continued until the cathode current is below 2 nA, which is in our experience the level 
of self sustained discharge. We have checked that in this way no sequentiality in the 
measurements is observed. 

The anode current is measured integrating over a time window depending on the pi- 
coammeter fuUscale. For the majority of the performed measurements the integration 
time of the picoammeter was ~600 /is. As a consequence the single current measurement 
is affected by large poissonian fiuctuations. The statistical error has been reduced to be 
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negligible, compared with the calibration and linearity error of the picoammeter, aver- 
aging over 300 single measurements of the current. Therefore we have assigned to each 
current measurement an error corresponding to 0.4% of the measured value, plus a zero 
point uncertainty of ±0.5 pA. Actually the LSD of the read out of the instrument corre- 
sponds to ±0.05 pA, but we have observed that the zero point current reproducibility is 
not at this level, likely due to some unavoidable parasitic current. 

From the values of the current we have obtained the gas gain dividing by the calculated 
primary current. Then we have corrected the calculated gain for the effect of small 
variations in the primary ionization current according to Eq. (jSj) using the pressure and 
temperature data, monitored for each measurement. From these values the estimated a 
coefficient has been obtained, by evaluating the derivative by discrete incremental ratios. 
In order to evaluate the robustness of the method against possible numerical instabilities, 
we have also checked that the discrete values obtained in this way are compatible with the 
smooth curve obtained from the differentiation of a polynomial fit to the data. Finally we 
have computed the effective anode reduced field and the value of the a coefficient using 
the actual pressure measurement, scaled to a reference temperature of 20° C. 

The entire cycle of measurements and the data reduction has been repeated six times, 
spanning from Ar-C02 100-0% to 0-100% in steps of 20%. We have observed that each 
set of curves is fully compatible inside the experimental error. In Fig. |21 which constitutes 
the central result of this work, we have reported the average of the values obtained in the 
five runs. 
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Pure gases 


Combined fit 


Ai 


(1.38 ±0.07) X 10^ 


(1.07 ±0.08) X 10^ 


Bi 


(6.05 ±0.08) X 10^ 


(5.6 ±0.1) X 10'* 


A2 


(2.2 ±0.9) X 10^ 


(2.6 ±0.2) X 10^ 


B2 


(2.6 ±0.3) X 10^ 


(2.65 ±0.05) X 10* 



Table 1 

In the first column the parameters to be inserted in Eq. [7| obtained from two separate fits 
of the data for pure gases only. The second column shows the results obtained from the 
fit of the model of Eq. © over the entire data set. Subscript 1 refers to CO2, and 2 to 
Ar. Errors are only the statistical error of the fit. 

4. Discussion 

Several parameterizations of the Townsend coefficient have been proposed in the liter- 
ature (for a review see e.g. Ref.p21 cind references therein). 
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Figure 3. Joint 95% confidence regions for tlie best fitted parameters relative to CO2. 
Tlie tliicker line refers to the results for the averaged data. 
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Figure 4. Joint 95% confidence regions for tlie best fitted parameters relative to Argon. 
The thicker hne refers to the results for the averaged data. 

However a simple thermal distribution of the electron energies, leads to a function of 
the type [12] 

a{S) = Ae-'^""^'^ (7) 

which has the appeal of a straightforward physical interpretation. In fact if we assume 
that the accelerated electrons have maxwellian velocity distribution with temperature 
kTe ~ eE Xcoii, the ionization rate from ground state of the gas is ^14] oc exp[—eion/kTe] 
giving B ^ (l/eXcoii), while A = {1/Xion) where Xcoii is the m.f.p. of the electron for 
elastic and inelastic scatterings, while Xion is the m.f.p. for ionizing scatterings only. In 
this context the () brackets indicate an average of the energy dependent m.f.p. over the 
relevant energy range of the electrons. 

It is remarkable that our data for pure Argon and pure CO2, reported in Fig. |21 can 
be well fitted with Eq. (0) in the lower part of the reduced field range, up to a value of 
a ~ 8, 000 ions/cm/atm. The fitted value are reported in the second column of Tab. H 
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Limiting our attention to mixtures of two gases, we propose here a simple additive 
model in which the ionization densities produced in the mixed gas is the sum of the 
densities in each gas. Therefore, starting from Eq. ^ for the pure gases, we can easily 
write for a mixture 

a=piAi e-^i^*/^ + p2 A2 (8) 

where B* = (pi-Bi +^2-82)5 Pi and p2 are the relative partial pressures, ei and 62 the first 
ionization potential of the two gases. 

In Fig. 121 we have reported as thin dot-dashed lines, the curves obtained from Eq. (jH)), in 
which we have inserted the coefficients A and B fitted on the data from pure gases only. 
It can be seen that the estimates of the Townsend coefficient for the various mixtures 
obtained in this way are not far from the measured values. 

More interesting is the result that we have obtained from a two-dimensional fit of the 
entire data set with the model of Eq. (jH)), leaving as free parameters the coefficients A 
and B of the pure gases, and considering the reduced field and the partial pressures 
independently measured coordinates of the data points. In Fig. El we show the result of 
the fit in this case as heavy dashed lines. It is clear that in this way we obtain a more 
reliable fit of the coefficients of the pure gases, because we use all the available information 
on the A and B coefficients for the pure gases at the same time. The fitted parameters 
are reported in the third column of Tab. Q We have also performed a separate fit to the 
curves obtained from of each run. As can be seen from Figs. El and |3J where the 95% 
confidence regions are reported, we also find that all of the fits of the single runs are 
compatible within 2a. 
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We conclude that our work shows that our proposed method of estimating the Townsend 
effective coefficient from the slope of the InM vs. InV^ curves is experimentally robust. 
It also shows that the simple functional form of Eq. (|7j) is adequate for predicting the 
evolution of the Townsend coefficient at high fields, if the detector is far from the regime 
of self-sustained regenerative discharge. In spite of its simplicity the additive model for the 
Townsend coefficient, in the form of Eq. (jH)) can be reliably used to predict the behavior 
of Ar-C02 mixture. It will be interesting to extend in the future this type of investigation 
to more complex mixtures, including gases with strong electron attachment and photon 
regeneration such as for example the CF4. 
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